Extreme times in financial markets 
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We apply the theory of continuous time random walks to study some aspects of the extreme value 
problem applied to financial time series. We focus our attention on extreme times, specifically the 
mean exit time and the mean first-passage time. We set the general equations for these extremes 
and evaluate the mean exit time for actual data. 
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I. INTRODUCTION 

The study of financial market microstructure has been 
the object of intense research since the beginning of 
mathematical finance [l|, Q ■ Thus one of the first mod- 
els of market microstructure was set by Bachelier who 
proposed the mechanism of the random walk (RW) to 
represent the dynamics of speculative prices Q]- In fact 
Bachelier's RW, with slight modifications and additions, 
has been successfully appHed to a variety of problems 
in finance ranging from price dynamics to option pric- 
ing 13, 3 ■ Unfortunately, RW models because of their 
Gaussian nature fail to account for a universal character- 
istic of price distributions: the fact that they are highly 
leptokurtic, with considerable fat tails at the wings of 
the distributions [a, Q . This was one of the main reasons 
which led Mandelbrot to propose the Levy distribution as 
an alternative description of the probability distribution 
of prices Q- However, pure Levy models have a serious 
drawback since any statistical moment, beyond the first 
one, does not exist. This fact has induced some authors 
to propose truncated Levy walks (TLW) as an alterna- 
tive model. TLW's have the advantage of having finite 
moments up to some order (which depends on the trun- 
cation procedure) while keeping, to some extend, some 
quite attractive properties of TL's such as self-similarity 
and scaling laws |3, E|- One of the problems of TL's is 
that the truncation is an ad-hoc procedure with little re- 
lation, if any, with the underlying dynamics. Alternative 
models based on the Gaussian process plus jumps [i3 or 
even continuous jumps 11] have also been proposed with 
acceptable results. 

Continuous Time Random Walks (CTRW's) are gen- 
eral models which perhaps better capture market mi- 
crostructure, specially that of high frequency data. 
CTRW's were first introduced by MontroU and Weiss in 
1965 il^] and have a long history of successful applica- 
tions to physics, chemistry and geophysics to name a few 
[ll[l3. To our knowledge, the application of the CTRW 
to finance is quite recent and it has not been fully devel- 
oped its potential [la,lla)ll3)ll3- One of the applications 
where CTRW's may represent a valuable achievement is 
in the field of risk control because the CTRW formalism 
provides a natural way of treating any particular real- 
ization of the price or return processes. This, in turn. 



facilitates the study of extreme events which are central 
to risk management. 

The statistics of extremes is a difficult field in prob- 
ability theory and its thorough description for a given 
random process can be quite involved, if not impossible, 
from an analytical point of view [l9ll2(il2ll |. In this paper 
we will use the CTRW formalism to study some aspects 
of the extreme value problem applied to finance. We will 
focus our attention on two of the simplest quantities re- 
lated to extreme statistics: the mean first-passage time 
(MFPT) and the mean exit time (MET). MFPT is the 
average time at which the random process reaches, for 
the first time, some pre-assigned value while MET is the 
mean time when the random process leaves, for the first 
time, a given interval. 

The paper is organized as follows. In Sec. II we outline 
the CTRW formalism and briefiy describe some of its 
applications. In Sec. Ill we develop the theory of MET 
for financial time series and apply it to real data in Sec. 
IV. Closing remarks are in Sec. V. 



II. OUTLINE OF THE CTRW 

In this section we summarize the main features of the 
CTRW formalism applied to the analysis of financial time 
series. We refer to the reader to [UO for a more com- 
plete account on the subject. 

Let S{t) be a financial price and let Iq be any 
initial time. The log-price or return is defined by 
Zit) = lnS'(i)/5'(io). If {Zit)) is the return mean value, 
we define the zero-mean return by 



X{t)^Z{t)-{Z{t)). 



(1) 



For the rest of the paper we will assume that the return 
is a stationary random process independent of any initial 
time io which is taken to be zero, io = 0- Hence t refers 
to a time interval. 

Suppose that X{t) is described by a CTRW. In this 
representation any trajectory consists of a series of step 
functions as shown in Fig. Q] Therefore, X{t) changes at 
random times t{),ti,t2,- ■ ■ ,tn,- ■ ■ and we assume that so- 
journs or waiting times, t„ — tn — in-i (^ = Ij 2, • • • , n), 
are independent and identically distributed random vari- 
ables described by a given probability density function 



(pdf) defined by ip{t)dt = Prob{t < t^ < t + dt}. 
At the conclusion of a given sojourn the return X{t) 
suffers a random jump described by the random vari- 
able AXn = X{tn) — X{tn^i) whosc pdf is defined by 
h{x)dx — Probjo; < AX„ < a; + dx}. We combine these 
two causes of randomness into one single density p(x, t) 
which represents the joint pdf of waiting times and ran- 
dom jumps, i.e., 

p{x, t)dxdt = Probja; < AX„ < a: + dx; t < Tn < t + dt}. 

We will further assume that p{x, t) is an even function of 
X so that there is no net drift in the evolution of X{t). 
Note that if waiting times and jumps are independent 
random quantities, then p{x, t) — h{x)'ip{t). In any other 
situation one has to specify a functional form of p{x, t) 
that is compatible with the observed data. Moreover, 
since the jump pdf, h{x), and the waiting-time pdf, ^{t), 
are the marginal densities of the joint density, any pro- 
posed form of p{x, t) must satisfy 



i,{t) 



p(x,t)dx; h{x) 



p{x,t)dt. (2) 



The main objective of the CTRW is obtaining the so- 
called propagator, that is, the probability density func- 
tion of the zero- mean return X[t): 

p{x, t)dx = Probja; < X(t) <x + dx}. (3) 

In [13, [iai we have obtained a general expression for the 
joint Fourier-Laplace transform of the propagator, 

p{uj, s)= e-'^dt / e-~''^''p{x, t)dx, 



in terms of the Laplace transform of the waiting time 
distribution, ^'{s), and the Fourier-Laplace transform of 
the joint distribution, p{uj, s). This expression reads 



p{uj,s) 



1 - piuj, s) ■ 



(4) 



Equation Q) constitutes the most general solution to 
the problem. In order to proceed further we would need 
to specify the form of the joint density p(a;, t) which has 
to be guessed from data 22] . There are nonetheless some 
general results that are independent of the model chosen 
for p{x,t) [13 • We briefiy outline them: 

(a) If the mean waiting time (r) is finite and the 
jump pdf h{x) has a finite second moment, 
P2 = (AX^) < 00, the asymptotic distribution of 
returns for long times approaches to the Gaussian 
density: 



p{uj,t) 



-M2^^t/2(r> (i>(r)). 



(5) 



and the asymptotic variance grows linearly with 
time: 



(X\t))^p^t (i»(r)). 



(6) 
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FIG. 1: A particular trajectory of the zero-mean return pro- 
cess along with a particular value of the random variable 

ta,b{xo). 



(b) lih{x) is a long-tailed density, i. e., h{x) ^ \x\ ^ " 
as \x\ — > 00. Then h{LLj) ~ 1 — fc|u;|" as w ^ for 
< a < 2. Moreover, if we assume that for ui 
small {re^^^-^) ~ (r), then the asymptotic return 
pdf approaches to the Levy distribution: 



p{uj,t) 



-,-fcli 



"*/<"> {f>{T)). 



(7) 



(c) At intermediate times, Ik, (r), the behavior of 
p{x,t) for large values of |a;|, is the same as that 
of the jump distribution: 



P{x,t) ^ -r^h{x). 
r) 



(8) 



From points (a) and (b) we see that the CTRW for- 
malism unifies two of the approaches mentioned in Sec. 
I: the Bachelier's RW (Gaussian models) and Mandel- 
brot's Levy models, thus providing a natural way of un- 
derstanding the truncation procedures appearing in TLW 
models. Finally, point (c) shows that fat tails in price dis- 
tributions have their origin in fat tails in the jump dis- 
tribution. Hence, by observing the jump pdf h{x) which 
is more easily accessible than the complete pdf p{x,t), 
one can infer the approximate behavior at the tails of 
p{x,t) via Eq. ^. We therefore see that the CTRW for- 
malism unifies in a single model the main characteristics 
observed in the empirical distributions of prices. 



III. EXTREME TIMES 

We assume that the zero-mean return process X{t) is 
described by a CTRW with a given waiting-time distri- 
bution, V'(Oi jump distribution, h{x)^ and joint density, 
p{x,t). Suppose that initially the return has a known 
value A'(O) = Xq [23. For t > we ask ourselves the 
following question: at which time interval X{t) leaves a 
given interval [a, b] for the first time? In other words, at 



which time the return is greater than certain value b or 
smaller than a for the first time? We call this quantity 
the exit time out of the interval [a, b] and denote it by 
ta.bi^o)- Obviously ta^b{xo) is a random variable since it 
depends on the particular trajectory of X{t) chosen (see 
Fig. ^. Our main objective here is to obtain, based 
on the CTRW formalism, the mean exit time (MET) 

The standard approach to MET problems requires the 
knowledge of the survival probability of the process in 
the interval [a,b] u^^^UM- Although the interest in 
knowing the survival probability is beyond any doubt, its 
attainment turns out to be quite involved. In this paper 
we present a direct and much simpler derivation of the 
MET and leave obtaining the survival probability for a 
later work. 

We decompose T{xo) in two summands 



T{xo) = Ti{xo) + T2{xo), 



(9) 



where ri(xo) is the MET to leave the interval [a,b] in 
only one jump and T2{xo) is the MET when more jumps 
have occurred. Note that in terms of the joint density 
p{x,t), Ti(xq) can be written as 

/■OC /"OC 

ri(a;o) ~ I tdt / p{x ~ xo,t)dx 
Jo Jb 



tdt 

vi-oo 



p{x — xq, t)dx, (10) 



where the first summand is the mean time for the random 
walker to escape through the upper boundary x = b, and 
the second summand is the mean time to escape through 
the lower boundary x — a. If the random walker has not 
exited the interval in the first jump the process will have 
attained at time t some value x € [a, 6] inside the interval 
and from that point the mean exit time will be exactly 
T{x). That is, 



T2(a;o) = / dt p{x - xo,t)[t + T{x)]dx. (11) 

Jo Ja 

Substituting Eqs. l(TT)|l - fTT|l into Eq. © and some sim- 
ple algebra involving the use of Eq. Q finally yield the 
following integral equation for the MET T{xo) 



T{xo) = (r) + / h{x - XQ)T{x)dx, (12) 

J a 

where (r) is the mean waiting time between jumps: 

(t) = / tiP{t)dt. (13) 



From a mathematical point of view Eq. (|12|l is a Fred- 
holm integral equations of second kind. Depending on 
the specific nature of the kernel h{x) there are some an- 
alytical approaches which allow to get an exact solution 



[25J . In the most general case if the kernel norm defined 

by 



h'^{x — y)dxdy, 



(14) 



is finite there is always a series solution that in many 
situations can be useful to obtain a good approximation 
[23. In the next section we will see examples of exact 
and approximate solutions. 

An important point should be emphasized: the fact 
that, as shown in Eq. H12|l . the MET does not depend on 
the possible coupling between waiting times and jumps. 
In other words, T(xo) is independent of the particular 
form of the joint density p{x, t). 

Another extreme time closely related to the MET is 
the mean first-passage time, MFPT, defined as the mean 
time at which the process attains a given value Xc for 
the first time. Let Ta^xo) be the MET out of the in- 
terval [a, b] then if xq < Xc the MFPT to Xc is defined 
as Tc(xo) — lima^_oo Ta^x^ixo), while if xq > Xc we have 
Tc{xq) = linife^oo T!r^.b(a;o)- Thus we see from Eq. (|12f) 
that, if Xq < Xc, the MFPT obeys the integral equation 

h{x - xo)Tc{x)dx, (15) 

-00 

with an analogous equation when xq > Xc- Unfortu- 
nately, Eq. H15(l is a singular integral equation for which 
the kernel norm defined above is infinite '2(?|. In such a 
case the solution of Eq. (|15|l may not exist ^]. Again, 
we will see an example of this in the next section. 



IV. SOME PROPERTIES AND RESULTS 

Equation H12|l constitutes the main result of this pa- 
per. In this section we will summarize some of the gen- 
eral properties of Eq. H12|) and present some results. We 
stress the fact, already mentioned in the previous sec- 
tion, that T{xo) depends on the entire jump distribu- 
tion, h(x), but does not depend on the whole waiting- 
time pdf, ipit), being only necessary for the existence of 
the MET that the mean waiting time (r) be finite. We 
also observe that Eq. H12|l is written for a general and 
non-symmetrical interval [a,b]. However, we can always 
transform the problem in order to work on a symmetrical 
interval, something that lightens the problem. Thus we 
can easily see that the symmetrical MET defined by 



tsym 



{x)=T{x + 



a + b 



(16) 



satisfies the equation 



rL/2 

rsym(a;o) = (t) + / h{x - Xo)Tsyra{x)dx, (17) 

J-L/2 

where L — b — a is the length of the interval. When 
the jump distribution, h{x) = h{—x), is even we see 



from Eq. |[T7|l that Tsym(— a;) satisfies the same equa- 
tion as Tsyjn(a;). In other words, Tsy^(x) = Tgymi—x) is 
also an even function which implies that the derivative 
of Tsym(2;) at X = is zero if it exists, T'-^{0) — 0. 

There are few cases in which it is possible to get an 
exact solution for the MET problem. This is the case 
when the jump pdf is given by the Laplace distribution 



h{x) = le 



(18) 



where 7 > and (AX^) = 2/7^ is the variance of jumps. 
One can show that in this case the integral equation H12() 
is equivalent to the following differential equation 



r'ixo) 



-l^r) 



(19) 



with boundary conditions 

T'{a)^^[T{a)-{T)l T' [b) = -^[T{b) - {r)]. (20) 



The solution to this problem reads 



T{x,) - Y 






7 2^0 



a + b 



(21) 

infinite since 

00. Conse- 



Note that in this case the MFPT is 
T{xo) — > cx), both as a — > — cxd and b 
quently Ti,{xo) — 00. We also note that the MET is a 
quadratic function of the interval length L = b — a. This 
is even more clearly seen by assuming a symmetrical in- 
terval b = —a = L/2 and that the initial return is zero, 
xo = (in fact, this is the usual situation [23 )■ In this 
case Eq. (I21II reads 



T(0) 



2 



(22) 



It is very illustrative to compare the above expressions 
for the MET with those of the ordinary random walk. 
If the price process follows a RW then in the continuous 
limit the zero-mean return is the Wiener process, i.e., 
X{t) = aW{t), where a is the volatility. In this case the 
MET out of [a, b] is E^ 



1 



Tnwixo) = -^{xq - a){b- xq). 
(J- 



(23) 



Observe that boundary conditions now are rRw(a) — 
Tuwib) = which are quite different than those give by 
Eq. (jSOJ. In order to compare this time with that of the 
CTRW just obtained in Eq. (|21|) we will scale both times. 
We thus define the following dimcnsionlcss MET's: 



and 



/TH* 



T^^ixo) = (aV2i2)TRw(a;o) 



/M = (1/(t)7 L )TcTRw(a;o) 



Then for a symmetrical interval b 
from Eqs. ^ and (|^ 



-'cTRwV'''OJ — -'rwv-^oJ ^" 2r2 o^T 



L/2 we get 
(24) 



where 



TnwM 



1 _ _Xo_ 

8 ~ 2L2- 



(25) 



We see that T*.^^^ 



{xq) -^ T*^{xo) when L > 7"^, that 
is, when the length of the interval is much larger than 
the jump volatility which is proportional to I/7. In this 
case our CTRW approaches the Wiener process with a 
volatility given by cr = \/2/j{t). 

Let us now obtain an approximate solution which will 
be valid for any sufficiently smooth kernel h{x). Suppose 
we have an even and zero-mean jump density satisfying 
the following scaling condition 



h{x) 






(26) 



where k > is the volatility of the jump process (note 
that the density (|18l) satisfies this condition with k = 
\/2/^). In this case Eq. (|17|l can be written as 



f{u) = 



L/2k 



L/2k 



H{v - u)T{v)dv, (27) 



where T(u) = Tsyui{Ku) and —L/2k < u < L/2k. Once 
we have an expression for T{u) the mean exit time is 
given, via Eq. ((T?)|l . by 



Tixo) = f 



2xo — a — b 
2k. 



(28) 



Suppose now that e = L/2k, is small. In this case an 
approximate solution to Eq. (|27|l can be easily obtained 
through an iteration procedure with the final result up 
to second order 27]: 

f{u) = (r) [1 -f 2H{Q)e + {H'{Q+) + 4H{0)^)e^ 
+ H'{0+y + O{e^)]. 



Hence 

nxo) = (r) 



l + H{0) 



[H'{0+)/4 + H{0f 



""'^'^xo-a-br + of^ 



Ak^ 



In the symmetrical case with xq — we have 



(29) 



r(o) 



l + H{0)[- 

K 



[H'i0+)/4 + H{0f] 



(30) 
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FIG. 2: Empirical distribution of the jump distribution h{x) 
for the U.S. doUar/Deutsche mark futures market. This figure 
suggests the presence of a power law which is adjusted by 
Eq. 13111 with the parameters reported in the main text. 



We see froni these expressions that, in general, the 
MET has for sufhciently small intervals a quadratic 
growth behavior as in the case of the exponential 
density H18|l . In fact, in this latter case k = \/2/7, 
H{x) = exp[-V2|a;|]/V2 and Eq. (jSHI) agrees with 
Eq. (|21|l . In other words the approximate expression be- 
comes the exact solution. 

Let us finally apply the above results to real data. 
The data consist of tick-by-tick prices for the U.S. dol- 
lar/Deutsche mark future exchange from January 1993 to 
December 1997 (a total of 1,048,590 data points). The 
empirical jump distribution is plotted in Fig. [3 We have 
shown elsewhere [l7l [l§ | that this empirical distribution 
is very well fitted by a power law of the form 



h{x) = 



(/3-1) 



27(1 + 1x1/7)'^' 



where (3 = 5.52 and 7 = 2.64 x 10""*. 
mean waiting time between jumps is (i 
rest of parameters are 

1V2 



(31) 

In this case the 
) = 23.65 s. The 



V(/3-2)(/3-3)' 



and 



HiO) 



(/3-1) 



V2(/3-2)(/3-3)' 



H'{0+) = 



/3(/3-l) 
(/3-2)(/3-3)' 



with numerical values k — 1.25 x 10 ^, H{0) = 1.07 and 
H'{0+) = -2.81. 

In Fig. Owe compare the empirical MET from the U.S. 
dollar/Deutsche mark data with the analytical approxi- 
mation given by Eq. 1)301) . We see there that when h{x) 
is estimated by the power law (|31|) . with the numerical 
values for k, H{0), and H'{0+) given above, the the- 
oretical expression does not fit the empirical data in a 




normalised semi-interval length L / 2k 

FIG. 3: Dots represent the empirical MET (measured in 
hours) as a function of the normalized semi-interval L/2k, 
fronr high frequency data of the U.S. dollar/Deutsche mark 
futures market. We also plot several theoretical predictions 
explained in the text. The best adjustment corresponds to Eq. 
13011 with all parameters directly estimated from data without 
making any hypothesis on the form of the jump distribution 
h{x). 



satisfactory manner. The reason for it lies in the fact 
the approximation for T{xo) given by Ea. (|30|l basically 
depends on the values of h{x) around x — while the 
power law density (|31() has been obtained to fit the tails 
of the jump distribution rather than its center. 

On the other hand, we can easily evaluate the values 
of K and H(0) directly from data without assuming any 
hypothesis on the form of the jump distribution. This 
yields K = 1.70 x lO"'', H{0) = 4.45 x lO'^. We also 
need an estimation for H'{0+) although this parameter 
is quite difficult to evaluate because data binning en- 
tails a certain degree of arbitrariness. Moreover the very 
existence of tick units questions the concept of deriva- 
tive. A first approach to solve this problem would be 
to evaluate H'{0+) using finite differences, which gives 
H'{0+) = 0.55. Nevertheless if we represent T(0) as 
given in Eq. (|30|l using the numerical values of k, H{0) 
and H'{0+) just reported, then the resulting function fits 
to some extend the empirical MET albeit being not com- 
pletely satisfactory. A better fit of T(0) can be obtained 
by using the same values of k and H{0) but estimating 
H'{0+) in order to furnish an optimal fit for the MET 
which finally yields the value H'{0+) — 1.54. This is rep- 
resented by the solid line curve of Fig. |3| Note that this 
final estimation shows that the approximate expressions 
for the MET given by Eqs. ^ and ^ are very effi- 
cient since they almost exactly reproduce the empirical 
behavior of the MET far away from the validness of the 
approximation itself, which we recall was only effective 
for small values of the normalized interval L/k. We can 
thus conclude that the MET presents a quadratic nature 
regardless the length of the interval. All of this seems 
to indicate that the approximations given by Eqs. I|29() 
and (|30|l have a more universal character than one might 
expect in advance. 

We finally observe that the Laplace density given by 



Eq. I^Hl: which yields the exact MET ^^, results in a 
good approximation of the MET for small values of the 
interval as we can also see in Fig. Olln this regard the 
Laplace density can be considered as a first approxima- 
tion to any more realistic jump distribution. 



V. CONCLUSIONS 

Using the CTRW framework for market microstruc- 
ture, we have developed a somewhat little known aspect 
of the problem: the study of extreme events, specially the 
mean exit time out of an interval, T{xo), and the mean 
first-passage time to some critical value, Tc{xo), although 
the latter turns out to be infinite in many situations. We 
have shown that these extreme times obey an integral 
equation which depends on the jump distribution h{x) 
and the mean waiting time (r). We have exactly solved 
the integral equation for the MET in the case where the 
jump distribution is governed by a Laplace (symmetric 
exponential) probability density function. We have com- 
pared this MET with that of the ordinary random walk 
model and showed that the MET for the Laplace density 
is bigger than Trw (the MET when X{t) is assumed to 
be the Wiener process). This seems to indicate that the 
models based on the Wiener process may underestimate 
the MET (specially for small intervals). In other words, 
RW models imply that the return process escapes faster 



from a give interval than models based on the CTRW. 
We believe that this can have practical consequences in 
risk control as well as in pricing exotic derivatives. 

We have also solved the integral equation for the MET 
using an approximate scheme which yields a solution 
valid for a general h{x) but when the length L of the 
interval is small (in an appropriate dimensionless scale). 
We have showed that for the Laplace form of h{x) the 
approximate solution for small L coincides with the ex- 
act solution, therefore showing that the Laplace is a 
"first-order approximation" to any sufficiently smooth 
and symmetrical jump distribution. We have applied 
the approximate solution to real data (high frequency 
data on U.S. dollar-Deutsche mark futures market) with 
a very good agreement between the empirical and the 
theoretical MET. 

We finally mention that more complete measures of 
the extreme events, such as the survival probability, with 
greater impact on risk control are under present investi- 
gation and some results will be published soon. 
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